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Abstract. We study the evolution of the spectral index and the normalization (flux) of the non-thermal component of the 
electron spectra observed by RHESSI during 24 solar hard X-ray flares. The quantitative evolution is confronted with the 
predictions of simple electron acceleration models featuring the soft-hard-soft behaviour. The comparison is general in scope 
and can be applied to different acceleration models, provided that they make predictions for the behavior of the spectral index 
as a function of the normalization. A simple stochastic acceleration model yields plausible best-fit model parameters for about 
77% of the 141 events consisting of rise and decay phases of individual hard X-ray peaks. However, it implies unphysically 
high electron acceleration rates and total energies for the others. Other simple acceleration models such as constant rate of 
accelerated electrons or constant input power have a similar failure rate. The peaks inconsistent with the simple acceleration 
models have smaller variations in the spectral index. The cases compatible with a simple stochastic model require typically a 
few times 10 36 electrons accelerated per second at a threshold energy of 18 keV in the rise phases and 24 keV in the decay 
phases of the flare peaks. 
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1. Introduction 

The intense hard X-ray emission observed during solar flares is 
the direct signature of the presence of highly energetic supra- 
thermal electrons. The quest for a viable particle acceleration 
mechanism drives both theoretical and observational investiga- 
tions. In his review on particle acceleration in impulsive so- 
lar flares Miller (1998 1 presents "the major observationally- 
derived requirements" for particle acceleration. The two main 
observational facts cited by Miller and used to put constraints 
on the electron acceleration mechanisms are the time scales of 
the acceleration (Is for acceleration from thermal energies to 
100 keV) and the acceleration rates (10 36 to 10 37 electrons s _I 
accelerated above 20 keV, to be sustained for several tens of 
seconds). Interestingly, the shape of the accelerated electron 
spectrum and its evolution in time are barely mentioned. Grigis 
& Benz (2004 henceforth Paper I) have analyzed X-ray obser- 
vations looking for systematic trends in the spectral evolution 
of 24 impulsive solar flares and confirmed the predominant 
soft-hard- soft (SHS) behavior of the observed photon spectra, 
first noted by Parks & Winckler ( 1969 1. This not only applies 
to the global evolution, but is even more pronounced in individ- 
ual peaks. The authors also give a simple quantitative descrip- 
tion of the SHS pattern, deriving an empirical relation between 
the normalization of the non-thermal component of the pho- 
ton spectrum and its spectral index. Can this systematic trend 
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in the evolution of the hard X-ray photon spectrum be used to 
put constraints on acceleration mechanisms? As pointed out al- 
ready in Paper I, the SHS behavior contradicts the idea that the 
flux evolves by a varying rate of identical, unresolved events, 
termed 'statistical flare' in avalanche models (Lu & Hamilton 
IT99TI 

Here we study the constraints the new quantitative infor- 
mation on the SHS behavior puts also on other acceleration 
models. Some conventional simple acceleration scenarios (pre- 
sented in Section are tested whether they can reproduce 
the observed spectral behavior and what constraints on their 
parameters can be obtained. Our goal is to demonstrate the 
method, to stimulate further comparisons between observation 
and theory and to call attention to the fact that the spectral evo- 
lution cannot be neglected by a successful acceleration theory. 
While the scenarios presented here are admittedly simple, we 
think that this first step will be extended in the near future to 
encompass more sophisticated models. 

The main piece of information that we use for this com- 
parison is the relation between the normalization of the non- 
thermal component of the spectrum (assumed to be a power- 
law) and its spectral index. Paper I studied this relation for pho- 
ton spectra. Here, we go a step further, and recover the electron 
spectra assuming an emission from a thick target, which also 
yields a power-law spectrum for the electrons. Therefore we 
can use the data from Paper I and convert the photon spectral 
index and normalization to the corresponding values for the 
electron spectra. This yields discrete time series of the spectral 
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index 5(t) and the power-law normalization O f0 (f) at energy £q. 
The acceleration models described in Section [5] provide the- 
oretical functions O f0 (5) depending on the model parameters, 
which can be fitted to the observed pairs (6(t), O f0 (f)). Applying 
this repeatedly for different flares and different emission peaks 
during flares, the distributions of the best-fit model parameters 
can be derived. 

We summarize in Section |2 the data reduction process 
yielding the dataset. It is used in Section|4]for the comparison 
with the models described in Section[3] 

2. Observations and Data Reduction 

We give here a brief summary of the data reduction process, de- 
scribed in full detail in Paper I. The photon spectral data used 
for this work is exactly the same as the dataset used in Paper I, 
where hard X-ray observations from the Reuven Ramaty High 
Energy Solar Spectroscopic Imager (RHESSI) of 24 solar flares 
of GOES class between Ml and XI have been analyzed. For 
each event count spectra were generated with a cadence of one 
RHESSI rotation (amounting to about 4 seconds). Each spec- 
trum was fitted with a model consisting of an isothermal emis- 
sion (bremsstrahlung continuum and atomic emission lines) 
characterized by a temperature T and an emission measure M, 
and a non-thermal component characterized by a power-law 
with index -y, normalization Fg a at energy Eq, and low-energy 
turnover at energy E lum . The full detector response matrix was 
used to calibrate the spectra. The fittings were done by means 
of an automatic routine, but were checked one by one, exclud- 
ing cases where the thermal and non-thermal emissions could 
not be reliably separated, or the non-thermal part was not well 
represented by a single power-law. After this selection we had 
a total of 91 1 good fittings for 24 events. For this work, we just 
use the time series of the spectral index y and the normaliza- 
tion of the power law F# for all the 24 events. The non-thermal 
component of the spectrum is thus approximated by 



in 



where F(E) is the photon flux at 1 AU in photons s _1 cm ~ 2 
keV _1 . We now go one step further than Paper I and transform 
the photon spectra into electron spectra. We choose the well- 
known analytically solvable thick target impact model (using 
the non-relativistic Bethe-Heitler cross section and collisional 
energy losses) to recover the injected electron spectrum. It is 
still a power-law 



(2) 



O(e) = O eo I - 



where O(e) is the total number of electrons s 1 keV 1 over the 
whole target. The electron spectral index is given by 



6 = y+l 

and the electron spectrum normalization is 
(6-1) (6- 2) 



<D £ =KF E Ele t 



o c o 



A<5-2,l/2) ' 



(3) 



(4) 



where f3{x, y) is the beta function, and the constant K is given 
by 

3n 2 e 4 In AD 2 



K 



6.4- 10 33 keVcm 2 



(5) 



(Brown [T971> . This transformation yields the time series <£> £o (f) 
and (5(f), which will be used for the comparison with the accel- 
eration models. 

The transformation of the observed photon spectrum into 
an electron spectrum can potentially alter the results of the 
comparison between observations and models, since it requires 
a knowledge of the physical conditions in the emission region, 
of their evolution in time and of the importance of the en- 
ergy loss processes. However, this is a necessary step for the 
comparison with acceleration models. For simplicity we have 
neglected the effects of photon reflection in the photosphere 
(albedo) and nonuniform target ionization. 

The assumption of a thick target leaves open the origin of 
temporal variations. In the following, we will assume that they 
are due to the acceleration rather than to a variable release of 
trapped particle. 

3. Acceleration scenarios 

We present here the frameworks of two acceleration scenarios, 
the constant productivity scenario and the stochastic acceler- 
ation scenario, which we will compare with the data on spec- 
tral flare evolution. A scenario can comprise different models. 
All models are required to yield a power-law distribution <1> 
in electron energy e (Eq. A specific model is defined by 
a unique functional relation between the spectral index 5 and 
the spectrum normalization <5 eo . The relation depends on some 
model parameters, which are assumed to be constant during a 
flare or a subpeak. This enables us to compare the model func- 
tion with the observed dataset (6(f), <E> f0 (f)). The first scenario 
is based on ad hoc assumptions, while the second is derived 
from a stochastic electron acceleration model considered e.g. 



by Benz (1977). Newer models along the stochastic acceler- 
ation line, like the transit-time damping model proposed by 



Miller at al. f 1996), do not imply a simple functional relation- 
ship, and therefore cannot be included at this stage. The spec- 
tral evolution of transit-time acceleration must be calculated 
numerically and is not yet available. Petrosian & Liu (2004i 
also give plots of electron spectra calculated from their model 
of stochastic acceleration by parallel waves, but do not pre- 
dict the time evolution of the spectra either. For a quantitative 
comparison between observations and theory, we are limited to 
models making concrete and usable predictions on the relation 
between spectral index and power-law normalization. 

3.1. The constant productivity scenario 

A class of models may be defined assuming that the productiv- 
ity of the accelerator is constant above a threshold energy e„. 
Either the electron acceleration rate N (electrons s~') to ener- 
gies above e» or the total power P input to electrons above e, 
are held constant, but the acceleration process evolves in such 
a way that flux and index vary. Although such processes have 
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Fig. 1. Electron energy spectra with spectral indices of 6 — 4, 6, 8 for each of the four models indicated on top and presented in 
Section|3] The vertical dashed line represents the lower energy threshold e, for the first 3 models, and the pivot point energy in 
the last plot. 



never been proposed, they may fit some of the data. Upon inte- 
grating Eq. (0, the models give the following relations between 
the electron spectral index 6 and the electron spectrum normal- 
ization & : 



N(6- 1) (e. 



— (constant rate model) 



P(6 - 2) / e, 



(constant power model) 



(6) 



(7) 



The total power in Eq. is expressed in the somewhat unusual 
unit of keV s _1 . The power in erg s _1 can be easily obtained 
multiplying P by the conversion factor 1.602 ■ 1(T 9 ergkev _1 . 

3.2. The stochastic acceleration scenario 

The index-flux relation in stochastic acceleration was explored 



by Benz (1977i and further elaborated by Brown & Loran 
(1985 1. In this model plasma waves accelerate stochastically 
the electrons in a plasma slab. From the diffusion equation for 



the electron distribution function Benz gets the approximate re- 
lation (corresponding to Eq. 20 in Benz. [T977l with the spectral 
index transformed from the electron number distribution into 
the corresponding value for the electron flux distribution) 



1 8e 2 In An V2e» 1 

6 = -- + + ■ =: -- + d + e , 

2 nW neyJW^L 2 



(8) 



where e is the electron charge, In A the Coulomb logarithm, n 
the ambient plasma electron density, L the length of the plasma 
sheet and Wo the time-dependant spectral energy density in 
the waves causing the acceleration. Following Brown & Loran 
( 1985 1, we recognize that the term d can be written as d = ae 2 
with 



47re 4 In AnL 



(9) 



where Lmfp is the Coulomb collisional mean free path of the 
electrons in the plasma sheet. The acceleration process needs 
Lmfp » L to be effective, and therefore a <£ 1. Benz ar- 
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Fig. 2. Plot of the spectral index 5 vs. the spectrum normalization <t>6o for different times in three flares (first column: 20 Feb. 2002 
09:54, second column: 04 Apr. 2002 15:29, third column 17 Apr. 2002 00:39). The points belonging to a certain rise phase are 
represented by stars, the points belonging to a decay phase are represented by triangles and the other points are marked by dots. 
The evolution of the constant rate and the stochastic acceleration model are represented by the dashed and full lines, respectively. 
In the upper row they were fitted to the rise phase points, in the lower row to the decay phase points. 



gues that the total electron flux <I>tot above e* is proportional 
to WqL, that is 

Otot = KW L. (10) 
Using Eqs. and d!0> we get the following relation <5 eo <-> 6 



*eA° C(0-1) 



6* 



1 + yjl +4a(d+ i) 



(ID 



where C = 2Ke 2 /(7i 2 e 2 ). In the collisionless case a — it sim- 
plifies to 



e.\ C 8-1 



2 ' 



(stochastic acceleration model) (12) 



Figure^shows the spectra for different values of 6 and their 
<D £o (5) according to Eqs. (|6j, (0 and (I12> . Interestingly, the 
model described by Eq. (I12> has the special property that all 
the spectra cross each other in a very narrow region of the plot, 
therefore exhibiting a behavior very similar to the one given by 
the pivot point model described in the following. 



The pivot point model assumes that all the non-thermal 
power-law spectra in the time series cross each other in an 
unique point. This pivot point is characterized by its energy 
e« and its flux <!>». The following relation holds between the 
spectral index 6 and the power-law normalization O f0 at energy 
£0 



*e,=<M- 

k«0 



(pivot point model) 



(13) 



The relation between 6 and i> f0 is given by Eqs. J(>J, 0, 
(II 1> . H3i for, respectively, the constant rate model, the con- 
stant power model, the stochastic acceleration model and the 
pivot point model. All models depend on two free parameters 
assumed to be constant during a flare or an emission peak. The 
threshold energy for the first 3 models and the energy of the 
pivot point in the last model have been all represented by the 
same symbol e», such that it is easier to compare the different 
equations. When more distinction is needed we will refer to 
them as ef ATE , e t POWEE \ e, ST0C and ef IV for the different mod- 
els. The second parameter, denoted N, P, C, O,, respectively, 
characterizes the flux normalization in the 4 models. 



P. C. Grigis and A. O. Benz: Spectral evolution of solar X-ray flares. II. 5 
CONSTANT RATE MODEL: RISE PHASES CONSTANT RATE MODEL: DECAY PHASES 




0.1 1.0 10.0 100.0 

e, [keV] 

STOCHASTIC MODEL: RISE PHASES 



10 44 - 

• 

10 42 - 
10 40 - 

,'(0, — 

O 10 38 - 
10 36 - 

10 34 - 

10 32 _ 
0.1 




1.0 10.0 
e, [keV] 



100.0 



10 44 - 

10 42 - .. 



1 10 38 




0.1 1.0 10.0 100.0 

£* [keV] 

STOCHASTIC MODEL: DECAY PHASES 



10 44 
10 42 
10 40 h 
O 10 ; 
10 

10 34 h 

10 : 
0.1 




1.0 10.0 
e. [keV] 



100.0 



Fig. 3. Best-fit model parameters for the stochastic acceleration and the constant rate model, separately for the rise and decay 
phases. Contour lines corresponding to 1<x,3<t (thick line) and 5<x levels for a Gaussian peak fitted to the data are also shown. 
The region under the dashed line has an electron acceleration rate lower than 10 37 electrons s _1 for events with spectral index 
<5 = 6. 



3.3. Relations with the pivot point model 

We now compute the energy of the approximate pivot point IV 
which results from a stochastic acceleration model given by 
the parameters (C, e^ T0C ) with a = 0. To compute the position 
of the pivot point, we find first the position of the intersection 
point of two spectra given by 5i,<Di = O^^i) and 62,^2 = 
•frfoOfe) m the stochastic acceleration model. Since the spectra 
are straight lines in logarithmic representation, using Eqs. (|2ji 
and H2i it is straightforward to find for the intersection e INT 

In &zl _ ? In S -^X 
eiNT ln<5 2 -ln<Di m fc-i s s + { 



In 



t STOC 



62 - S\ 



62 - 8\ 



(14) 



To find the approximate pivot point, we take the limit of the 
previous expression for 62 — > S\ . Putting 6\ — 6 and 62 = 5 + A 
with the condition that A/ 6 <SC 1, we get to the first order in A/ 6 



In 



EiNT 
~STOC 



5- 



(tf-l)(tf+i) 



(15) 



Therefore two spectra with spectral index around 5 will have a 
common point at 



_ JTOC 
e INT - £, 



exp 



(<y-D(tf+i)J 



,STOC 



f(S). 



(16) 



The function f(8) depends weakly on delta for 5 > 3, and there- 
fore all the spectra will cross each other in a narrow energy and 
flux range. The energy of the corresponding pivot point is ap- 
proximately ef IV = ef T a e t ST0C • /(6) = 0.90 e» ST0C . Similar 
relations holds approximatively also for the other models: 

1 



e PIV s ef ATE -exp 



5- 1 



epiv 



e, POWER ■ exp 



1 



6-2 



(17) 



(18) 



In the constant productivity models, the pivot point has an 
energy which is slightly larger than e„, and in the stochastic 
acceleration model the pivot point energy is slightly lower than 
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e». Therefore in these models most accelerated electrons have 
energies comparable with the pivot point energy. On the other 
hand, the purely phenomenological pivot point model does not 
request the presence of electrons at energies close to the pivot 
point energy: the pivot point may be virtual in the sense that 
the electron energy distribution may be a power-law at higher 
energies and turnover at an energy e T 



^TURN > f PIV 



4. Comparison with the data 

4.1. Fits in 5-® eo plane 

We proceed now to compare the observed evolution of the 
spectra with the models presented in the previous section. The 
dataset described in Section[2]consists of discrete time series of 
the electron spectral index <5(f,) and of the power-law normal- 
ization O f0 (f,) at energy eo at the times t\, t%, . , . ,t n for the 24 
different flares. For each model described in Section[3]we have 
a relation between O f0 and 5, and the comparison of the data is 
done by a least-square fitting of the inverse function <5(€> £o ,P) 
to the observed pairs (<D £o (f,), (5( ?,-)), where P is the model pa- 
rameter vector. Since the spectral index 6 varies over a smaller 
factor than the flux normalization cD eo , it is better to use the lat- 
ter as independent variable for the fitting, and therefore we use 
the inverse function <5(<£> eo ) instead of 0)^(6). The choice of eo 
is arbitrary, but we settled for eo = 60 keV, consistent with the 
observed range of electron energies. Paper I remarked that the 
scatter of the data is smaller on average for time series belong- 
ing to rise or decay phases of the non-thermal emission peaks. 
A given rise or decay phase is defined as a series of at least 
three consecutive points showing an increase or a decrease of 
the flux. Therefore we independently fit the model parameters 
to each rise and decay phase. Figure |2 shows (^>(,q,6) points 
consequent in time measured for three flares. The constant rate 
and stochastic acceleration models corresponding to the best fit 
to the rise and decay phases are both shown in all graphs. In 
each fit the model parameters are constant in time. The curves 
corresponding to the two different models are nearly identi- 
cal because Eqs. l|6} and dl2> are functionally similar. The two 
other models (constant power and pivot point) yield curves (not 
shown in Fig.|2ji which also are much alike the ones shown. 



4.2. Comparing model parameters 

As a next step we compare the best-fit parameters between 
the different models. Is one of the models preferable, yielding 
more plausible values? Note that there is an evident correla- 
tion between the two parameters of each model. It is due to 
the fact that the points (6, <&6o) of each phase lie approximately 
on straight lines with different slopes, but all of them passing 
near the point (<t>6o = 10 31 electrons s _I kev _1 ,<J = 6.5). This 
geometrical constraint in the 6-<&(,q plane is the cause of the 
correlation between the two model parameters: if e, is small, 
N, P, C, or <t>„, respectively, must be large. The distributions 
of the best-fit model parameters for the constant rate and the 
stochastic acceleration model are presented in Fig.[5]separately 
for the rise and decay phases. Since the parameters vary over 



several orders of magnitude, we chose a logarithmic represen- 
tation. 

In all the four models, we found that most of the parame- 
ters are concentrated in a relatively narrow region, except for 
about 20-30% of the points lying far away from the peak of 
the distribution. The average value and standard deviation of 
the distribution are poor estimators of the position and width of 
the central peak in such a case. Instead, we characterized the 
distributions by fitting two-dimensional Gaussian functions to 
the data. This procedure yields a reasonable estimate for the 
position and width of the distribution's peak. To further ensure 
that outliers do not have strong influence, we restricted the fit 
to the region where e, > 0. 1 keV and the second fit parameter 
is smaller than 10 45 (in each parameter's units). The number of 
cases in the region to be fitted is given in Tabled and its range 
is displayed in Fig. [3] The model parameter distribution must 
be binned for the fitting procedure. We chose bins centered on 
the median values of the two parameters and with a width of 
0.5 times the average deviation from the median. 

In Fig. |3 the contour lines corresponding to the lcr, 3cr 
(thick line) and 5cr levels of the two-dimensional Gaussian dis- 
tribution are superimposed on the data. We define the outliers 
as being the points outside the 3<x contour line. The number 
of outliers is about 23 % and does not vary significantly in the 
different models and phases (Table0. 

The outliers correspond to rise or decay phases where the 
spectral index changes little as the flux increases or decreases. 
In the framework of the pivot point model: if the pivot point 
for the spectrum lies at very low energy and high flux, it is 
possible to have large variations at high energies with small 
changes in the spectral index. The other 3 models behave in a 
fashion similar to the pivot point model (as shown in Fig.Q, so 
this explication also applies for them. 

With the exception of the pivot point model, the parame- 
ters of the other models have a direct physical meaning. We 
can check if they lie in an acceptable range when compared 
with other observations and generally accepted values on solar 
flares. For the sake of simplicity and uniformity, we derive first 
the total power injected in the accelerated electrons (in erg s~') 
and the acceleration rate of electrons (in electrons s -1 ) for the 
different models. The following equations express the electron 
O(e) de as a function of the model 

parameters: 



POWER . 

P(6 - 2) 
e*(<5-l)' 

C 



N 



STOC 



(£+i) 2 ' 



(19) 
(20) 

(21) 



The following equations express the power P — k e <I>(e) de 
as a function of the model parameters: 



pRATE = KN£t 



6-1 
6-2' 



pPOWER _ K p 



(22) 
(23) 
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Table 1. Results of model fitting to the observed evolution of the spectral index and non-thermal flux in the rise and decay phases 
of the non-thermal emission peaks. 



Model 


Equation for 


Phase 


Nr. of 


Nr. of phases 


Nr. of phases 


Fitted 


Centroid of the 


Half width of 




Oa, «-> 5 




phases 


in fit region 


in 3-tx region 


variable 


distribution 3 


the distribution 13 


Constant 


o 


Rise 


70 


64 


52 


N 


5.70 • 10 35 el. s" 1 


15 


rate 












€* 


13.9 keV 


1.6 






Decay 


71 


68 


56 


N 


2.05 • 10 35 el. s _I 


8.4 














6* 


18.6 keV 


1.5 


Constant 
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Rise 


70 


64 


51 


P 


1.18 • 10 3 ' keV s- 1 


8.3 


power 
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13.6 keV 


1.6 






Decay 
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67 


55 
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5.85 ■ 10 36 keV s _1 


4.5 
















17.6 keV 


1.4 


Stochastic 
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c 


8.09 • 10 36 s" 1 


16 


acceleration 














18.1 keV 
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56 


C 


2.23 • 10 36 s" 1 
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24.2 keV 


1.4 
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53 


3>, 


9.39 ■ 10 34 el. s" 1 keV 
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e, 


17.5 keV 
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Decay 


71 


67 


56 


<D : . 


1.28 • 10 34 el. s _1 keV 


1 11 














e 


22.3 keV 


1.4 



(a) centroid of the best fit 2-dimensional Gaussian distribution 

(b) expressed as a multiplicative factor 
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Fig. 4. Distribution of the electron acceleration rate (left column) and the total power in the accelerated electrons (right column) 
for the constant rate model (upper row) and the stochastic model (lower row). The distributions are shown in gray for the rise 
phases and black for the decay phases. The tail of the distribution to the right of the main peak contains events with unphysically 
high values for the rates and powers. 
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->STOC 



6-1 



(5-2){5+\) 2 " 



(24) 



where k = 1.602 ■ 10~ 9 erg kev _1 . Figure |4] shows the distri- 
bution of the acceleration rate and the electron power for the 
constant rate model and the stochastic model, for both the rise 
phases (gray) and the decay phases (black). The distributions 
show a central peak and an extended tail to the right. Rise and 
decay phases in the tail represent about the same subset in all 
models, requiring extremely high values for both acceleration 
rate and power. This is due to the fact that these phases fit mod- 
els having a very low threshold energy e», and the integral of 
the spectrum diverges in the limit e» — > 0. The low value of e* 
is necessary to account for the events where the spectral index 
varies slowly during the phase. 

What are the maximum values acceptable in reality? The 
interpretation of hard X-ray observations of solar flares re- 
quires a large number of accelerated electrons. In his review 
on the flare mechanism, Sweet (1969 1 already requires 10 36 
electrons s to account for the observed hard X-ray emission. 
Brown & Melrose t ! 1977b derive a requirement of 5 ■ 10 36 elec- 
trons s _I accelerated above 25 keV. Miller ( 119981 cites 10 36 - 
10 37 electron s above 20 keV. Recent RHESSI observations 
of X class solar flares yield similar values (e.g. Holman et al. 
2003 Saint-Hilaire & Benz 2005). If the non-thermal electron 
spectrum extends to energies lower than 20 keV, this number 
could be higher. This usually does not contradict the obser- 
vations, since during a large flare, the photon flux at energies 
lower than 20 keV is typically dominated by thermal emission. 
From a theoretical point of view the number of flare electrons 
and energy available in the active region environment is lim- 
ited. Some mechanism for electron replenishment needs to be 
operative but is often left unspecified. Therefore it is not clear 
how the sustainable electron acceleration rate is limited. As 
an upper limit for comparison with our results, we take as the 
highest reported electron acceleration rate N MAX = 10 37 elec- 
tron s . There is also a physical limit on the power that can be 
injected into the accelerated electrons. As a generally accept- 
able maximum power p MAX , we will assume 10 29 erg s in 
agreement with the above cited authors. 

Some of the derived model parameters are highly implausi- 
ble, requiring acceleration rates and power input far exceeding 
the above limits. In Fig. [3] the limit for 6 = 6 is indicated in 
the N-e* and C-e* planes. Table Ogives the number of cases 
in total, as well as in Fig. [3] and in the region used for fitting 
a Gaussian distribution. Figure [5] indicates that most rise and 
decay phases below the limit of plausibility are within the 3-<x 
limit of the distribution, and thus corroborate the Gaussian fit. 
Its centroid and half-widths (in log presentation) are given in 
Tabled 

Since we follow the time evolution of the spectra from the 
onset of the non-thermal emission, there is a large difference 
of more than 2 order of magnitude between the value of the 
minimum and maximum normalization <l>6o. Surprisingly, the 
clustering of phases in the 7V-e, and C-e* planes (Fig. [3j indi- 
cate that most flare peaks can be interpreted by models within 
a relatively small range of threshold energies e«. The second 
model parameters N, P, C, O, are spread upon a larger range. 



They include the effect of the different sizes of the non-thermal 
emission peaks of our sample. Note that e* is smaller in the rise 
phase than in the decay of a flare peak on average. On the other 
hand, the normalization is larger in the rise phase on average. 

Gan ( 1999) reported the presence of a pivot point in the 
SMM/GRS spectra of two X-class flares. His results of 41 keV 
and 77 keV for the energy of the pivot point of the electron 
spectra are larger than ours. This may be due to the fact that he 
analyzed larger events, or that he fitted energies mostly above 
our energy range, thus possibly above an high-energy break in 
the spectrum. 

The value of the spectral wave energy density Wo in the 
stochastic acceleration model can be computed. From Eqs. 
(0, (ED and C = 2Ke 2 /(n 2 e 2 ) we get: 



W L = 



2el 



n 2 e 2 (6 + i) 2 



(25) 



Using the average values of e, = 20 keV and 6 = 6.4, this 
amounts to 



W 2 x 10" 



10 7 cm 



erg cm 



(26) 



This value is somewhat lower than the ones cited in Benz 
(1977) and Brown & Loran (1985), because we have used a 
lower value of e* and a larger value of 6. In the stochastic 
model, the total energy density Uj = J W(k) die in the turbulent 
waves can be computed by 



U T - — ^ 3 x 10 



10 1() cm 



— erg cm 



(27) 



where Ad is the Debye length. It is well below the magnetic 
energy density, which is larger than unity. 

5. Discussion 

The four models presented in Section [3] are based on different 
assumptions, but have similar functional relationships between 
the electron spectral index 6 and the non-thermal spectrum nor- 
malization O f0 , as required by observations. As a consequence, 
it is possible to visualize the spectral evolution of all these mod- 
els by the presence of a common pivot point in the non-thermal 
electron spectra during the rise or decay phase of each hard X- 
ray peak. This general property however, inhibits selecting the 
best model from spectral fits. 

The comparison of the observed evolution of the non- 
thermal spectrum during rise and decay phases of a non- 
thermal emission peak with the models yields best-fit model 
parameters. Considering the simplicity of the assumed mod- 
els it is surprising that the distribution of the model parame- 
ters is reasonable for about 77% of the 141 observed rise/decay 
phases, but yields an unphysically high electron acceleration 
rate for the rest of the observed events. The numbers of outliers 
is not statistically different between the models. 

The events in the unphysical region of parameter space 
manifest themselves as the enhanced tail of the two- 
dimensional distribution shown in Fig. [3] The simple models 
do not succeed in describing these events, which all exhibit a 
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slowly varying spectral index. These models would totally fail 
to reproduce events with constant spectral index. The mathe- 
matical reason of the failure is that these model just assume 
power-law behavior of the observed electron distribution above 
a fixed threshold energy e„. To account for the events with 
slowly varying spectral index, the fit converges toward very 
low values of e„, whereas the models can be justified only in 
the case where e, is larger than about a few keVs. More sophis- 
ticated models, like the one proposed by Miller et al. (0996 » 
or Petrosian & Liu (2004i, explicitely accelerate non-thermal 
electrons out of a thermal distribution. They do not show a 
power-law behavior at low energy and therefore should be less 
susceptible to such low-energy divergencies. 

The best-fit model parameters show an asymmetry between 
the rise and decay phases of the emission peaks. Such an asym- 
metry was already reported in Paper I. In the constant produc- 
tivity scenario, this would correspond to a reduced productivity 
of the accelerator in the decay phase. Similarly, the fitting to the 
stochastic acceleration model suggests a lower electron accel- 
eration rate on average for the decay phase. 

6. Conclusion 

We compared observations (described in Paper I) from 
RHESSI on the evolution of the normalization and spectral in- 
dex of the non-thermal component during emission peaks of so- 
lar flares with model predictions by means of fittings the model 
parameters to the observed values for the electron spectral in- 
dex 5 and the spectrum normalization (Deo at 60 keV (as shown 
in Fig. |3. All the models selected for this study, described in 
Section 13 feature the soft-hard-soft behavior for a single flare 
peak and fit well the observed spectral evolution of flare peaks. 
We have shown in Section [3] that they have a functional de- 
pendence between 5 and <t>6o which closely resembles the one 
given by the pivot point model (Ea.ll3>. For this reason the data 
cannot discriminate between models. 

While it is possible to fit reasonably well the models to 
the observed data, the resulting model parameters imply un- 
physically high electron acceleration rates and energies for the 
20-30 % events were the spectral index changes slowly as the 
flux rises, thus showing less prominent soft-hard-soft behavior 
(Fig. 0}. This is due to the fact that the models need a very low 
threshold energy e, (that is, the energy above which electrons 
are accelerated, corresponding approximately to the energy of 
the pivot point) to account for low rates of change of the spec- 
tral index. A low value for e* enables the model to provide rea- 
sonably good fits to the data, but will often lie outside the range 
of validity of the models studied. For instance, the stochastic 
acceleration model will break down for electron energies com- 
parable to the thermal energy of the medium in which the elec- 
tron are accelerated, since the interaction of the accelerating 
waves with the thermal component was not considered for this 
simple model. The assumption of a fixed e, during a peak may 
also be challenged. Allowing e» to change during the observed 
events could reduce the excess acceleration rate values. 

The quantitative relation between the observed spectral in- 
dex and the normalization of the spectrum exploited here has 
the advantage that it does not depend explicitely on the time 



evolution of the two variables. Therefore it is comparatively 
easier to compare with model predictions than the full-fledged 
time evolution of the spectrum. The 5 <-> <1> relation puts an 
additional requirement on the acceleration mechanism. 

The comparison between observations and models pro- 
posed here can be done with real data in a straightforward way, 
and it does provide interesting results for relatively simple the- 
oretical models: a majority of flare peaks (about 77 % of the 
rise and decay phases) can be well fitted by these simple mod- 
els within a compact region of parameter space. Nevertheless, 
the rest (outliers with unphysical model parameters) point to 
the fact that more degrees of freedom are necessary for inter- 
pretation. While there is a wealth of data available thanks to the 
RHESSI mission, we lack concrete predictions by more com- 
plex acceleration models on the behavior of the spectral index 
during flares. 
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